Neutron scattering evidence for isolated spin-1/2 ladders in (C5Di2N)2CuBr4 
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Inelastic neutron scattering was used to determine the spin Hamiltonian for the singlet ground 
state system of fully deuterated BPCB, (C5Di2N)2CuBr4. A 2-leg spin-1/2 ladder model, with J± = 
(1.09 =b 0.01) meV and Jy = (0.296 =b 0.005) meV, accurately describes the data. The experimental 
limit on the effective inter- ladder exchange constant is \Jint\ ^ 0.006 meV, and the limit on total 
diagonal, intra-ladder exchange is |Jf + Jf' \ < 0.1 meV. Including the effects of copper to bromide 
covalent spin transfer on the magnetic form-factor, the experimental ratios of intra-ladder bond 
energies are consistent with the predictions of continuous unitary transformation. 

PACS numbers: 75.10.Jm, 75.40. Gb, 75.50.Ee, 78.70.Nx 



I. INTRODUCTION 

Low dimensional quantum magnets, with integer spin 
per unit cell, frequently exhibit a macroscopic singlet 
ground state. While impervious to weak fields, there is 
a critical field beyond which an extended critical phase 
can exist, and it is of considerable interest to explore the 
spin dynamics at the quantum critical point and within 
the putative critical phase. 

The n-leg antiferromagnetic (AFM) spin ladder con- 
sists of n parallel chains of magnetic moments with AFM 
exchange interactions along the chains (Jy) and between 
neighboring chains (J_l). In keeping with the Lieb- 
Schultz-Mattis theorem,!^ the excitatio n sp ectrum has 
a gap for n even and is gapless for n odd ! ^*^ * ^ * The n = 2, 
S = 1/2 case is of particular interest because such ladders 
may form dynamically in copper oxide superconductors 
and play a role in suppressing magnetism in favor of su- 
perconductivity. The Hamiltonian for the simple two-leg 
spin-ladder is given by the first two terms in 

H = J\\ ^ S_^-,^ • S_^-+i,/ + J± ^ S_^-,i • S_^-,2 

j,l • Sj_i,2 

3 3 

-\-J ^ ^ ^j,m ' ^j,n J ^ ^ Sj,m ' Sji + l,n 

j,'m,n j,'m,n 

-g^sH S,,, , (1) 

where / = 1,2 indexes each of the two chains and j is 
the rung index. However, for BPCB (see Fig. [ija)) ad- 
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ditional interactions are possible. The third and fourth 
terms are frustrating diagonal intra-ladder exchange in- 
teractions, Jf and J^?/, between spins in neighboring 
chains and rungs. The next two terms describe two pos- 
sible inter-ladder interactions, J' and J'', where m and 
n denote adjacent chains in different ladders. We define 
= J' — J" to discuss the effective interaction between 
adjacent ladders. The last term in the Hamiltonian is the 
Zeeman term associated with an applied magnetic field, 
H. It is included here in anticipation of future high field 
experiments, but for this work, H = 0. 

Consider the ideal ID system (J', J" = 0). If Jy = 
and one of the frustrating exchanges is zero, H describes 
an alternating spin chain, with the physics controlled by 
/3 = Jf/J±- Another extreme, Jf = Jf' = 0, is the 
ideal spin ladder, where the physics is controlled by a = 
J\\/J^. In the limit |a| oo, the system is composed 
of decoupled ID chains with a gapless spectrum.^ ^ ^ 
Any finite produces an isolated singlet ground state 
with a gap A ^ \Ji.\/2 (see Ref. 5). This general state 
of affairs persists into the strong coupling limit, |a| <C 1, 
where the ground state is a singlet, separated from the 
lowest lying triple t of excited states by an energy gap 

Given the extensive theoretical and numerical treat- 
ments that Eq. ([l]) has received,^ it is of great interest 
to identify materials containing spin-1/2 ladders that can 
be driven to quantum criticality through the application 
of a magnetic field. Excluding the system described in 
this paper, spin ladder materials known thus far either 
have energy scales that are too large to be affected by 
an applied field^^ or have significant inter-ladder inter- 
actions that i nduce Neel order above the critical field.^^ 
DiagonaP^^ and cyclic^^ exchanges, along wit h the pos- 
sible coexistence with other magnetic systems,'^^^^^EII20] 
represent additional challenges that are of interest in 
their own right but that may complicate analysis of the 
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quantum critical phenomena. 

Here we report experimental evidence that fully 
deuterated BPCB, (C5Di2N)2CuBr4, is a nearly ideal 
realization of an assembly of non-interacting spin- 1/2 
ladders. We provide an accurate determination of the 
relevant exchange constants, and quantify the low tem- 
perature exchange bond energies, in the absence of any 
applied magnetic field. These values can be used for de- 
scribing the field induced quantum critical state in other 
experiments. 



II. EXPERIMENTAL DETAILS 

We used inelastic neutron scattering to 
probe magnetic excitations in fully deuterated 
bis(piperidinium)tetrabromocuprate(II) , commonly 
referred to as BPCB, (C5Di2N)2CuBr4. BPCB is 
monoclinic (space group P2i/c), with room temperature 
lattice parameters a = 8.49 A, 6 = 17.22 A, c = 12.38 A, 
and P = 99.3°.^ Throughout we shall denote wave 
vector transfer in the corresponding reciprocal lattice 
c\{hkl) = hsi* + /cb* + /c*. The Cu^+ ions form ladders, 
as shown in Fig. [l] The legs of the ladders run along a, 
and the rungs are nearly along c*, with a small tilt of 
24° above the a - c plane.^^^J^ 

High field magnetization measurements {{) < H < 
30 T) were performed by Watson et al.^^'^ revealing a 
lower critical field of Hd = 6.6 T, an upper critical 
field of Hc2 = 14.6 T, and an infiection point at half 
the saturation magnetization. Through careful compar- 
ison of bulk thermo-magnetic data to various models, 
BPCB was identified as a two- leg spin ladder in the 
strong coupling limit with J^/ks = 13.3 K (1.15 meV) 
and J\\/kB = 3.8 K (0.33 meV). The analysis indicated 
that BPCB possesses an isolated singlet ground state for 
H < Hci and forms a gapless Luttinger spin liquid for 
Hci < H < Hc2' 

Inelastic neutron scattering is a sensitive probe of 
atomic scale correlations and interactions in singlet 
ground state systems. Using this technique, we find that 
BPCB is highly one-dimensional {\J^^^/J± \ ^ 5 x 10~^), 
making it an excellent candidate for future exploration 
of the high field critical phase. The findings extend, and 
are consistent with, previous experimental results on the 
spin Hamiltonian for BPCB. 

Using 99.9% deuterated starting materials, the sam- 
ple used for our neutron scattering measurements was 
made by the same process and team of scientists as de- 
scribed in Ref. 22, The sample consists of five deuterated 
single crystals, with a total mass of 3 g and co-aligned 
within one degree, for scattering in the {hOl) reciprocal 
lattice plane. The measurements were performed at the 
NIST Center for Neutron Research using the time-of- 
fiight Disk Chopper Spectrometer (DCS).^^ The chop- 
per cascade was phased to provide an incident wave- 
length A = 5 A and energy resolution AE ~ 0.1 meV. 
The sample was cooled in a liquid helium cryostat to 




FIG. 1: (Color online) Various projections of the previously 
determined structure of BPCB.^^ (a) Exchange interactions 
entering the Hamiltonian in Eq. (IT]). Black lines indicate the 
legs and rungs of the ladder, (b-d) BPCB crystal structure, 
except hydrogen, projected along each of the crystallographic 
axes. 



T = (1.4 ±0.1) K. The non-magnetic background was 
measured at T = 25 K, where magnetic scattering is 
widely distributed in energy and momentum. To a good 
approximation, the high temperature scattering can be 
treated as independent of sample orientation. This back- 
ground measurement was subtracted from the 1.4 K 
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FIG. 2: (Color online) False color image of raw = 5 A 
neutron scattering data measured with c*||ki. Intensity is in- 
tegrated in the / direction, (a) T = 1.4 K. (b) T = 25 K, used 
as background, (c) Intensity at T = 1.4 K after background 
subtraction. For presentation purposes only, th e data in (c) 
were smoothed as described in the text (see Sec. IIIB). 



data.f^ The DAVE software package was used to perform 
the initial analysis, to compute the energy resolution, and 
to extract the data required for advanced processing.!^ 
An example of raw data measured with the c* axis 
parallel to is shown in Fig. [2|a). An integration over 
the k and / directions is performed to generate this fig- 
ure. Such a procedure is commonly used for analyzing 
data acquired on low dimensional systems using chopper 
spectrometers.^^ However, we will use our measurements 
to quantify inter-ladder exchange. For this purpose, all 
components of momentum transfer in the horizontal scat- 
tering plane were used in the subsequent analysis. The 
trajectories in the (/i, 0,/) plane are shown in Fig. [s] for 
energy transfers of 0.8 meV and 1.5 meV, and for all 
sample orientations used (^ = Zk^, c* = 0° , —10°, and 
+60.9°). The trajectories sample the (/i,0,/) plane suffi- 
ciently to test for the presence of inter-ladder dispersion. 




FIG. 3: (Color online) Trajectories in (/i, 0, /) plane for the 
three sample orientations measured at huo = 0.8 meV and 
1.5 meV, with = 5 A. 



III. ANALYSIS AND DISCUSSION 
A. Exchange paths 

As a first step towards a model spin Hamiltonian, 
we shall discuss the structure and chemical bonding in 
BPCB. A review of the magnetic exchange interactions 
in a wide-range of tetrabromocuprate^^ provided guid- 
ance on possible exchange paths. While these consid- 
erations are not rigorous, they can provide a reference 
against whic h to compare the experimental results. It 
was prop osecP^'^ that the rung interaction {Ji_) is asso- 
ciated with overlap of Br~ adjacent to copper sites, while 
the exchange interactions along the legs of the ladder 
(J||) are mediated by a combination of hydrogen bonds 
and non-overlapping Br~ orbit als, and therefore should 
be weaker. However, Fig. [T] shows that the Br-Br dis- 
tances associated with inter- and intra- chain interaction 
are in fact quite similar. Furthermore, these distances are 
4 — 5 A, which is more than twice the covalent radius of 
bromine. This observation indicates that magnetic inter- 
actions in BPCB are mediated by intervening hydrogen, 
as becomes increasingly clear when piperidinium radicals 
are included in the picture, Fig.[TJb-d). Excess hydrogen 
is located around the nitrogen sites in the piperidinium 
rings, and these are found aligned with the nearest ap- 
proach of bromine atoms, associated with neighboring 
Cu^+ ions. Two piperidinium groups are involved in pro- 
ducing Jx, and only one for Jy. The greater number of 
Br-H bridges for rung over leg interactions leads to an 
expectation of dominant rung exchange. 

Any frustrating diagonal interaction {Jp or Jf>) would 
involve traversing the piperidinium molecule. Note that 
a large \Jf\ > Jy and Jpf = would result in an al- 
ternating chain, as opposed to a ladder. This situation 
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is found, for example, in MCCL!^ Furthermore, Jp and 
Jp' being associated with different bond lengths of 8.96 A 
and 12.64 A, respectively, suggests Jp dominates. 

The strongest inter-ladder interaction, J', is expected 
between ladders separated by c, and it is mediated by hy- 
drogen bonding through the same piperidinium molecule 
involved in the leg exchange. In addition, a J" inter- 
ladder exchange interaction is possible between atoms in 
ladders separated by c + a. In conjunction with J', a 
finite J" interaction might produce frustration, giving a 
small J^, and reduce inter-ladder dispersion. Any inter- 
ladder exchange in the b direction would involve a longer 
path, through two piperidinium molecules, and is there- 
fore expected to be weak. 



B. Global Fitting based on Single Mode 
Approximation 

The single mode approximatioiP^ generally provides 
an excellent description of the dynamic spin correlation 
function for gapped quantum magnets. The assumption 
that all spectral weight resides in a dispersive "triplon" , 
combined with the first moment sum-rule, leads to the 
following expression for <S(q, hw): 
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FIG. 4: (color online) Bandwidth over bandgap (VK/A) ratio 
for spin excitations versus normaliz ed in tra-ladder exchange 
interactions, computed using Eqs. (4j5). The experimental 
result of VK/A = 0.73 for BPCB constrains intra-ladder ex- 
change interactions to lie on the dashed red line. 



5(q,M = -^ ^^^^ ^-^^ ^e,[l-cos(q.d)] ,(2) 

where is the triplon dispersion relation and 6d = 
Jd(So • Sd) are the so-called bond energies, which sum 
to the ground state energy for T = 0. The summation 
is over spin pairs with finite exchange interactions. For 
the Hamiltonian in Eq. ([T]), in each unit cell there is one 
J_L rung spin-pair term, two Jy terms, one Jp term, one 
inter-ladder J' term, and one inter-ladder J" term, as 
shown in Fig.[TJa). The dispersion relation, £^q, is a func- 
tion of J^, OL — J\\/J^^ P = Jp/J±, and 7 = 
The neutron scattering intensity is obtained by multi- 
plying 5(q, huj) by the square of the Cu^+ magnetic form 
factoi^^ and convoluting the result with the instrumental 
resolution. 

The raw data binning results in an effective wave- 
vector resolution approximatively equal to the pixel size. 
The energy resolution was wider than the pixel size, and 
it is described by a Gaussian, with the width depend- 
ing on the configuration of the instrument and on the 
incident and scattered neutron energies. All subsequent 
analysis focuses on scattering data in the 0.7 meV to 
1.7 meV energy range. Global fits to the subtracted 25 K 
background were performed with all three sample ori- 
entations simultaneously. Any residual background was 
treated as momentum and energy independent. To gen- 
erate false color images of the treated data for Fig. |2] 
(and Figs, [g] and [7| which will appear subsequently), 
a 3D Gaussian smoothing was used with full-widths at 
half-maximum of 0.03 meV and 0.03 r.l.u. along the h 
and / directions. 



1. Triplon dispersion 

To quantify the magnetic interactions in BPCB, 
Eq. JTl), we examine the triplon dispersion relation. Fig- 
ure [2[suggests that a lowest order approximation to the 
dispersion is a periodic lattice sinusoidal function of the 
form 

W 

£;q = A + y (1 + cos(27r/i)) + A^* cos(27r0 , (3) 

where A denotes the gap in the absence of inter-ladder 
dispersion, W is the intra-ladder bandwidth, and Ac^ 
is the amplitude of dispersion along c*, resulting from 
^int- found that this last term is zero, within experi- 
mental uncertainty, but a more detailed analysis will be 
employed later to establish an experimental limit. The 
global fit (Fig. [^f)) yields a bandwidth W = (0.62 ± 
0.03) meV, and a spin gap A = (0.85 ± 0.01) meV. The 
latter value can be compared with A = 0.82 meV, ob- 
tained from magnetization measurementJ^on a hydroge- 
nous powder sample, where Hd = 6.6 T and (g) = 2.13. 
The structure factor indicates that the s± bond energy 
is dominant. The next step towards a spin Hamiltonian 
for BPCB is to relate the phenomenological parameters 
characterizing the dispersion relation to exchange con- 
stants. 

For a strictly one-dimensional model {J'^J" = 0), 
we can use perturbative expressions for the dispersion 
relation to extract exchange constants from the data. 
Contributions to dispersion from Jp and Jpf cannot be 
distinguished when both are small, so we define Jp = 
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{Jf + JfO/2- When both a = J\\/J± and p = Jf/J± 
are present, the model is a ladder with frustrating diag- 
onal exchange, or equivalent ly, an alternating spin chain 
with next nearest neighbor exchange. In either case, the 
dispersion is given hy^^ 



— \^ a^(a, cos(27rm/i) 



(4) 



m=0 



with 

ao 

ai 



a 3 , 



= a 



2 ^ A 



{a- 1)3 



2^ 



(5) 



1/ /^^2/-, /^N 



To account for inter-ladder exchange, in the first ap- 
proximation, we add 7COs(27r/) to Eq. Q. For stronger 
coupling, or different inter-ladder exchange paths, one 
expects more complicated / dependence, and possible 
cross terms involving both h and / dependence. In the 
absence of J" exchange, = J' / J^. If J" is present and 
frustrates J', 7 = = {J' — J") / J^- In princi- 

ple there should also be a q-independent term associated 
with Jf^. This effect is, however, sufficiently small to be 
neglected. 

With Jx dominant, there are two limiting cases: J\\ — 
or Ji? = 0. When Jy =0, the system is an alternating 
spin chain, with dispersion controlled by = JfI Ji_^^ 
From the a\ term, we note that BPCB can be de- 
scribed only by a negative so Jf would have to be 



25 



20 



15 



10 









15 

P 

^ 14 
13 


















-0.1 0.0 0.1. 
P 



-1.0 



-0.5 



0.0 

p=Jp/J, 



0.5 



1.0 



FIG. 5: ifc2 as a function of Jf/J±, for T4^/A = 0.73. The 
dashed lines correspond to Hci — 13.8 T and Rci — 14.6 T, 
as reported in different references.^^^^ Several other measure- 
ments of the upper critical field lie in this range.^^ ^ ' ^^ ' ^^ ' 
Detail shown in the inset indicates that \ JfIJi\ < 0.05, so 
IJfAsI < 0.4 K. 
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FIG. 6: (Color online) Global fit to a spin ladder model, 
(a) color map of smoothed (see text), background-subtracted 
data, with = Zk^, c* = —10°. The color scale is the same 
one used in Fig. 2c; (b) color map of the corresponding fit; 
(c)-(h) constant energy cuts, with an energy window of 0.1 
meV. The solid lines represent the global fit to the data when 
including inter-ladder coupling effects. 



ferromagnetic, corresponding to a FM/AFM alternating 
spin chain. However, calculations show that one cannot 
achieve the experimental bandwidth over bandgap ratio, 
W/A = 0.73, in the small /3 limit where Eq. ([5| is valid. 
In addition, a strong ferromagnetic interaction is incom- 
patible with magnetization measurements. 

A contour map of the W/A ratio as a function of a 
and P is shown in Fig. [4j The dashed line corresponds to 
all (a, P) pairs that are consistent with W/A = 0.73. For 
all points on this line, with \P\ < 0.5, ao varies by less 
than 3%. Hence, any set of (a, /3) on this line accurately 
describes the observed dispersion relation, with the value 
of J± within 3% of 1.09 meV. An additional constraint is 
therefore required to uniquely determine (a, P). This is 
provided by previous measurements of the upper critical 
field ff^^, 2 2 32 . 33 , 34 , 35 36 ^ gp-^ ladder with frustrating 
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FIG. 7: (color online) False color images of neutron scattering intensity for three sample orientations (rows): (column a) 
smoothed (see text), background-subtracted, data; (columns b-e) simulations with the phenomenological dispersion relation, 
Eq. ([s]), but only one bond energy (see top) modulating the intensity for each column; (column f) fit including all bond energies. 
For each sample orientation, select sections of data with unique values of h and / are shown. The color scale is the same one 
used in Fig. 2c. 



diagonal exchange, the upper critical field is 

gi^BHc2 = J± + 2J|| , (6) 

where Hc2 is independent of J^, if Jf < Ji. and Jp < 
2J||."^Fi gure 5 shows Hc2 versus Jf/J^ for values (a, 
(3) along the dashed line in Fig. |4] The values found 
experimentally vary between 13.8 T (Ref. [S^) and 14.6 
T (Ref. "22), as indicated by dashed lines in the figure 
and in the insert. From this analysis, we conclude that 
\Jf/Ji.\ < 0.05 or \JF/kB\ < 0.4 K. As an alternative 
method, we included the values of the upper critical field 
as a constraint in the fit to neutron scattering data and 
this yields /3 = Jf/ J± = —0.02 ± 0.10. These results are 
consistent with the upper limit reported by Watson et 
al.^ based on magnetization measurements performed 
at 700 mK. On the basis of this tight limit, the frustrat- 
ing exchange and the corresponding bond energy were 
neglected in the subsequent analysis of scattering data. 

When Jf is neglected, Eq. ^ describes the dispersion 
for an ideal ladder. A global fit to the scattering data 
yields Ji_ = (1.09±0.01) meV, Jy = (0.296±0.005) meV, 
and 7 = (0.002 ± 0.006). The quoted error bars reflect 
systematic error estimated as 10% of the energy resolu- 
tion. The statistical errors reported by the fitting routine 
were a factor of 2-3 smaller. 



The fitted q and E dependent intensity calculated 
for this model is shown in Fig. |6|b). A portion of 
the data, together with the fit, is presented as sev- 
eral constant energy cuts in Fig. [6]^c-e). The values for 
J± and J|| are in excellent agreement with the values 
obtained from neutron scattering j ^^ l ^^ l magnetization,!^ 
NMR,^ magnetostriction,'32ISI and specific heat and 
mametocaloric efFecPl measurements. We note that 
data has been fitted to an expression valid to third 
order of a and /3, Eqs. Q and ([5|. In principle, 
much higher order expressions can be obtained using 
other methods, including a particle conserving contin- 
uous unitarv transformation (QUTW^^^ and linked- 

I ' I I ^ 

cluster-expansion methods.^^1^^ Given the small val- 
ues of a and including higher order terms in the fitted 
dispersion is unnecessary for BPCB. 

2. Exchange bond energies 

Within the single mode approximation, Eq. ([2|, the 
exchange bond energies, = Jd(So • Sd), modulate the 
neutron scattering intensity periodically in q-d. A simu- 
lation of the intensity pattern, considering each possible 
exchange path in isolation, is shown in Fig.[7| where each 
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row corresponds to a different sample orientation. The 
observed h dependence of the intensity (column a) closely 
resembles that associated with the J± bond (column b). 
The J|| term yields a periodic modulation of intensity 
that is dissimilar to the data, and all other terms have 
intensity minima where the data has maxima. 

The fit to Eq. ([2| (see results in Table [ij) yields an un- 
reasonably large bond energy for the inter-ladder dimer 
considering the weak interactions. A possible explana- 
tion is that it is not appropriate to use an isotropic spin- 
only form factor for the Cu^+ ion.^^ Hubbard and Mar- 
shall showed that covalent bonds strongly affect neutron 
scattering intensities and magnetic form factors .^^ For 
the particular case of CuBr." anion, EPRfound,E^and 
calculations confirmed, that the electron density is sig- 
nificantly shifted from the copper d-orbitals into the a 
ligand orbit als. In the absence of a calculated magnetic 
form factor for BPCB similar to the one for cuprate spin 
chains,!^ we chose to modify the Cu^+ ionic form factor 
by isotropically rescaling its q dependence by a factor r 
to account for the spin density transfer to bromine. 

The global fit (Table shown in Fig. W\ (column f), 
yields e\\/ei_ = (0.05 ± 0.02), Ef/si, = (0.02 ± 0.03), 
and sj'^jnje^ = (±0.07 ± 0.15). A fit using the bond 
energies of both diagonal exchanges is not a significant 
improvement compared to these results, and the bond 
energy ratios presented above are essentially unchanged. 

Theoretical bond energies were obtained using a 
particl e conser ving continuous unitary transformation 
(CUT). 141142143] rj.^^ CUT is realized perturbatively at 
the isolated rung dimer limit. The elementary excita- 
tions conserved after the transformation are triplons.^^ 
The static correlation functions for perpendicular, par- 
allel, and diagonal bonds of the two-leg ladder can then 
be determined from the ground-state energy per bond, 
Eo{J±, J\\^ JF)/N\)onds^ which we calculated exactly up 
to order 7 in J\\/Ji_ and Jp/J^- Note that we are using 
the bare series in the following since we restrict the dis- 
cussion to small and intermediate values of the couplings 



TABLE I: Fit results with an ionic magnetic form factor and 
with a modified form factor that accounts for covalency ef- 
fects. For comparison, we show results from other neutron 
scattering experiments 





Ionic 


Covalent 


Other results 




form factor 


form factor 




Ji_ (meV) 


1.09 ±0.01 


1.09 ±0.01 


1.13±0.0l'^ 


J\\IJ^ 


0.288 ±0.057 


0.272 ± 0.002 


0.252±0.03^ 


Jf/J± 


-0.02 ±0.10 


-0.02 ±0.10 






-0.003 ± 0.003 


0.002 ± 0.006 


0.007^ 




0.10 ±0.02 


0.05 ±0.02 






-0.02 ±0.04 


0.02 ±0.03 




£j',j" /e± 


-0.27 ±0.04 


±0.07 ±0.15 




r 


1.00± 0.00 


2.36±0.13 






1.303 


1.293 
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FIG. 8: (color online) Solid black lines show bond energy 
ratios for Jf and J±, as a function oi [3 — Jf / for several 
values of a = J\\/ from 0.05 to 0.5, every 0.05, as estimated 
by the continuous unitary transformation. The experimental 
value and corresponding errors are shown as horizontal lines. 
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FIG. 9: Ratio of leg to rung bond energies versus a — J\\/ J±j 
as calculated using the continuous unitary transformation. 
Experimental values are shown as horizontal/vertical solid 
lines, with error-bars indicated by dashed lines. 



J|| and Jp. Using the Feynman-Hellman theorem, one 
finds for the static dimer correlation function 



C± - (S_^-,iS_^-,2) 



1 



bonds 



d_ 

±dJi 



in) 



(7) 
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with analogous expressions for other bonds. 

For comparison to the experimental data, calculations 
of Sp/^i. versus [3 = Jp/Ji. are presented for various 
values of a = J\\/ in Fig. [sl (black solid lines). The ex- 
perimental result Sp/^i. is shown as a red line in Fig. [sj 
with the dashed lines indicating uncertainty. This mea- 
surement does not impose additional constraints on Jp- 

For an ideal ladder, without frustrating or inter-ladder 
interactions, the leg to rung bond energy ratio versus 
a, computed by CUTs, is shown in Fig. [9j The bond 
energy ratio extracted from the neutron scattering data 
is shown with dashed lines. Given the rough nature of 
our approximation to the covalent form factor, the level 
of agreement is acceptable. 
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1.35 

1.30 



1.25 



-0.2 



-0.1 



0.0 



0.1 



0.2 



FIG. 10: Reduced for the global fit to BCPB inelastic 
neutron scattering data as a function of 7 = Jt^t/J^- The 
dashed line is a quadratic fit for points around the minimum, 
which indicates that 7 = 0.002 ± 0.006. 



3. Inter-ladder exchange 



The fit used in the previous section finds that J^^ is 
at least two orders of magnitude smaller than Jy and J^. 
To evaluate the robustness of this finding and obtain a 
quantitative uncertainty limit on 7, the fit was repeated 
for several different fixed values of 7. The resulting val- 



ues for the reduced are plotted versus 7 in Fig. [10 
A quadratic fit close to 7 = 0, yields 7 = 0.002 ± 0.006, 
and this analysis indicates at least two orders of mag- 
nitude difference between and Jy. Note that this 
result is specific to the assumed nature of inter-ladder 
dispersion. 



A second, less model dependent approach, involves 
generating constant q-cuts (width 0.05 r.l.u.) through 
the experimental data at different h values for all three 
sample orientations. The difference between the Gaus- 
sian fitted peak position to such cuts and the strictly 
one-dimensional dispersion relation is plotted as a func- 
tion of / in Fig. 11 The error-bars in the figure are the 



positional uncertainties of these Gaussian fits and should 
be compared to the energy resolution of the instrument, 
~ 0.1 meV. For some cuts, no signal was observed 
above background, in which case no point is shown in 
Fig. 11 There is no apparent systematic deviation from 
a zero residual as a function of again indicating the 
absence of magnon dispersion perpendicular to a*. 

Low temperature NMR measurements by Klanjsek et 
alW reveal 3D magnetic order for T < 100 mK. From this 
observation, an average inter-ladder coupling of ~ 20 mK 
(~ 1.7 /ieV) was inferred, and an identical result was 
found using neutron diffraction at high magnetic fields. 
These results considered four nearest neighbors in their 
mean field expansion. Therefore, the strength of the total 
effective inter-ladder exchange energy is ~ 80 mK. This 
value is comparable to the limits on | ^ 70 mK set 
by our fit. 




0.0 0.5 
1 (r.l.u.) 



FIG. 11: The difference between the best fit one dimensional 
dispersion relation and gaussian fits to cuts through raw data 
resolved versus wave vector transfer along c*. No systematic 
/ dependence is observed indicating that to within error a one 
dimensional model is appropriate for BPCB (7 = 0.002 ± 
0.006). 



IV. CONCLUSIONS 

We have shown that (C5Di2N)2CuBr4 is an excellent 
realization of two- leg spin- 1/2 ladder in the strong cou- 
pling limit. The inferred rung exchange Jx = (1.09 ± 
0.01) meV and leg exchange J\\ = (0.296±0.005) meV are 
in excellent agreement with values obtained from other 
techniques. Using two different methods of anal- 
ysis, we showed that the effective inter-ladder exchange 
is more than two orders of magnitude smaller than 
J^. These results confirm the previous conclusion that 
BPCB can be classified as a one-dimensional system. 
Alone, the neutron data do not provide a direct mea- 
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surement of Jp- However, in combination with high field 
magnetization studies, NMR measurements, and theo- 
retical calculations of the dispersion relation, the neu- 
tron data sets an upper limit on Jp + Jf'^ which is 
one order of magnitude smaller than the rung exchange 

{\{Jf + Jf')/Ji.\<Q.I)- 

The single mode approximation provides an excellent 
account of the data and the intra-ladder bond energies 
extracted are in agreement with results from continuous 
unitary transformations within experimental error. The 
intensity pattern can be understood only if covalency ef- 
fects are taken into account. The experiment elucidates 
the spin interactions in BPCB for analysis of recent^ 
and future high field neutron scattering experiments. 

During preparation of this manuscript, we became 
aware of the field dependent measurements of Thiele- 
mann et alW^ Our data at zero field confirm their re- 
sults, adding information regarding frustrating interac- 
tions and bond energies. 
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